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Spin up in RX J0806+15 - the shortest period binary 
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ABSTRACT 

RX J0806+15 has recently been identified as the binary system with the shortest 
known orbital period. We present a series of observations of RX J0806+15 including 
new optical observations taken one month apart. Using these observations and archival 
data we find that the period of this system is decreasing over time. Our measurements 
imply / = 6.11 x 10 -16 Hz/s, which is in agreement with a rate expected from the 
gravitational radiation for two white dwarfs orbiting at a given period. However, a 
smaller value of / = 3.14 x 10~ 16 Hz/s cannot be ruled out. Our result supports the 
idea that the 321.5 s period is the orbital period and that the system is the shortest 
period binary known so far and that it is one of the strongest sources of constant 
gravitational radiation in the sky. Furthermore, the decrease of the period strongly 
favours the unipolar inductor (or electric star) model rather than the accretion models. 
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1 INTRODUCTION 

For accreting stellar binary systems in which the secondary 
star (mass donating) is on the main sequence, the minimum 
orbital period is ~80 mins (Paczynski & Sienkiewicz 1981). 
For secondary stars not on the main sequence (such as he- 
lium stars), systems can also have orbital periods less than 
80 min (Savonije, de Kool & Van den Heuvel 1986 and refer- 
ences therein). The other option for short period systems is 
that the secondary star must be degenerate forming a dou- 
ble degenerate binary with a compact primary. A number of 
these systems have been known for some time although it is 
only relatively recently that their nature has been revealed 
(cf Warner 1995). The orbital periods of these systems lie in 
the tens of mins range. 

Recently a group of three binaries have been discov- 
ered in which the orbital period is thought to be less than 
~10mins: ES Cet (KUV01584-0939) 10.3 min (Warner & 
Woudt 2002), V407 Vul (RX J1914+24) 9.5 mins (Cropper 
et al 1998) and RX J0806+15, 5.4min (Ramsay, Hakala & 
Cropper 2002a, Israel et al. 2002). If these periods have been 
correctly identified as the binary orbital period, then these 
have the shortest known binary periods. As such they should 
be strong sources of gravitational waves. 



In the case of ES Cet, strong He lines are seen in optical 
spectra (Wegner et al. 1987) which have recently been found 
to be double peaked (P. Woudt, priv. comm.) suggesting 
the presence of a disk. In the case of V407 Vul, no emission 
lines are seen (Ramsay et al. 2002b) while RX J0806+15 
shows only weak He lines (Israel et al. 2002). Very recently 
(Israel et al. 2003) demonstrated that the optical and X-ray 
emission of RX J0806+15 are in anti-phase, as in the case 
of V407 Vul (Ramsay et al. 2000). 

Currently there are several models to explain the obser- 
vations of V407 Vul and RX J0806+15. A neutron star pri- 
mary and an intermediate polar model cannot be excluded 
although they are considered unlikely. The most likely mod- 
els all have a white dwarf primary and secondary. Two 
models (a polar and a direct accretor model) are driven by 
mass transfer while one (a unipolar inductor or electric star 
model) is driven by an electrical current. Of these models 
the unipolar inductor model fits the observational properties 
better than the accretion models, although these cannot be 
ruled out (Cropper et al. 2003 (review), Wu et al. 2002, 
Marsh & Steeghs 2002, Ramsay et al. 2002a and Norton et 
al. 2002). 

The best method to differentiate between these mod- 
els is to measure the period change of these systems. For 
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systems driven by mass transfer the period should increase 
over time (eg Savonije, de Kool & Van den Heuvel 1986, 
Nelemans et al. 2001). In the case of the unipolar inductor 
model, the system should be driven entirely by gravitational 
radiation and the period should decrease (Wu et al. 2002). 

Strohmayer (2002) showed using ROSAT observations 
that for V407 Vul the period is decreasing at a rate consis- 
tent with that predicted from the emission of gravitational 
radiation. In this paper we present observations of the 5 
min system RX J0806+15 and determine the change in the 
period. We then discuss the implications of this result. 



2 OBSERVATIONS 

In order to measure possible period changes in RXJ0806+15, 
we have obtained data from the Nordic Optical Telescope La 
Palma between 5th and 8th of January 2003 and 5th and 7th 
of February 2003. These observations consist of high speed 
CCD photometry with an approximate time resolution of 
15 sec (10 sec integration time). The observations were car- 
ried out using ALFOSC in imaging mode without any filter. 
We binned the CCD chip by a factor of two in both direc- 
tions and selected only a small 100x100 subwindow for fast 
readout. A total of 4687 images were obtained. These were 
bias corrected and fiatfielded in the usual manner. Condi- 
tions in January 2003 were photometric and the seeing was 
typically better than 1". In February, however, seeing was 
poorer (about 1".5). 

Additionally some archival ESO VLT data from 2001 
November 12 was analysed to provide further constraints for 
a period from our earlier (January 2002) NOT run. There 
were a total of 214 data points that were combined with our 
1314 datapoints from two nights in January 2002. These 
data were also reduced using standard procedures. All of 
the data were heliocentric corrected. 

There are also some archival Chandra observations 
available. However, given the length of the observation (20 
ksec), we have not included it here. 
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Figure 1. The Lomb-Scarglc power spectrum of the Jan-Feb 2003 
data and (inset) the same plot zoomed in near the best period. 
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Figure 2. The fits to the change in period assuming two different 
values for the ROSAT period. 



3 PERIOD ANALYSIS 

We have three different datasets that can be used to search 
for the change in the orbital period (Table 1). Chronologi- 
cally these are: I) The ROSAT observations from 1994-95, 
II) The combined VLT data from Nov 2001 and NOT data 
from Jan 2002 and finally III) The NOT data from Jan-Feb 
2003. We will now proceed to estimate the best period for 
each of these datasets and then use the results to derive the 
change in the orbital period. 

The earliest of our three datasets is the 1994-1995 
ROSAT data (Burwitz & Reinsch 2001). They found two 
possible values for the best period in their ROSAT data, 
namely: 321.5393 or 321.5465 s (± 0.0004 s). Their pref- 
erence was the first of these. For consistency we have re- 
analysed the ROSAT data and confirm their results. 

For our optical datasets we have used the Lomb-Scargle 
power spectrum (Scargle 1982) to search for the periods. 
However, in order to obtain reliable error estimates (crucial 
for the P measurement) we have performed Monte Carlo 
simulations of our time series. This enables us to take into 



account the true effects of aliasing in a rigorous manner. 
Our Monte Carlo approach consists of following steps, i) 
Find the best period using the Lomb-Scargle algorithm, ii) 
Fix the period and fit the light curve shape with a 2nd or- 
der Fourier fit. iii) Then use the (Gaussian) errors from the 
Fourier fit to produce 10000 synthetic datasets using the 
same time points as in the original dataset and add noise 
comparable to the noise in real data, iv) Finally analyze the 
10000 datasets around the best period. To do this in a rea- 
sonable amount of CPU time, we have used the minimum 
string length method (Dwortesky, 1983) which can be eas- 
ily restricted to the immediate neighbourhood of the best 
period. 

Our second dataset, which consists of VLT data from 
Nov 2001 and NOT data from Jan 2002 is problematic in 
terms of aliasing. Both of these runs are too short to 'phase' 
them together in a unique manner and as a result severe 
aliasing arises . However, we can use some a priori informa- 
tion from our later, Jan-Feb 2003, NOT run to remove the 
aliasing (i.e. pick the right spike from a handful of possi- 
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Observation 



Dates (HJD-2400000.0) 



Best period (s) Error (s) 



ROSAT (1994 &1995) 

VLT (Nov 2001) & NOT (Jan 2002) 

NOT (Jan -Feb 2003) 



49648.9-49658.2 

49822.9-49827.0 

52226.258-.373 

52289.533-.653 

52290.354-.740 

52645.443-.674 

52646.453-.639 

52647.435-. 773 

52676.420-.438 

52677.493-.576 



321.5393 OR 321.5465 



321.53352 



321.52957 



0.00040 
0.00034 

0.00045 



Table 1. The times of observation and best periods for our three different datasets 





Figure 3. The January 2002 (top) and January-February 2003 
white light light curves folded on the orbital period. 



ble spikes suggested by the Lomb-Scargle analysis and the 
subsequent Monte Carlo simulation). 

The P from Jan 2003 data combined with the ROSAT 
runs is about 1.2 msec/year. Further, our best period from 
NOT 2003 run is 321.52957s ± 0.00045s (see below). As- 
suming these values we can estimate that the period in 
Jan 2002 should be approximately 321.532s. The possible 
spikes present in the periodogram are 321.515s, 321.534s and 
321.552s. Only one of these (321.534s) has a period close to 
the ROSAT and NOT 2003 periods. If we now assume this is 
the 'correct' power peak, we can use our Monte Carlo study 
to get 321.53352s ± 0.00034s for our second run. 

The third run consists of data from three consecutive 
nights in Jan 2003 (5th-8th) and from two nights in Feb 
2003 (5th and 6th). This run does not suffer from serious 
aliasing and using the Monte Carlo technique outline above, 
we find 321.52957 ± 0.00045 s as the best period. We have 
plotted the periodogram of the NOT 2003 run in Figure 1. 
The inset shows the part near the best period. 



4 THE PERIOD CHANGE 

To measure the possible change in the orbital period. We 
have fitted our periods from three different epochs (defined 
as the mean of the time points from each of the epochs) with 



a linear fit (Figure 2). Depending on the choice of which 
ROSAT period we use, we obtain two sets of frequency (/) 
and / (= -P/P 2 ). 

If we adopt 321.5393s as the true ROSAT period (as pre- 
ferred by Burwitz & Reinsch, 2001). We get / = 3.14 x 10~ 16 
Hz/s. On the other hand, the longer ROSAT period 321.5465 
s yields / = 6.11 x 10~ 16 Hz/s. In both cases the error is 
2.0 x 10~ 17 Hz/s. Judging from our fit, the latter period of 
the two is better. Neither of the fits is very good, but given 
the number of measurements to be fitted (three) , we cannot 
make a definitive conclusion between the two solutions. 

We have also folded both the 2002 and 2003 NOT 
data over the orbital period and plotted that in Figure 3. 
In order to do this, we need to adopt a value for / that 
gives us the correct phasing. This can be used as a fea- 
sibility test for our model where the period is decreasing. 
We find two epochs when the flux should be at minimum, 
one just before the 2002 NOT dataset and another before 
the 2003 NOT dataset. These are HJD: 2452288.999256 and 
2452645.436696. Now, to produce correct phasing an integer 
number of (changing) periods must have passed in between 
these two epochs. As we expect the period change to be lin- 
ear, at least in short time scale, we can estimate an average 
number of periods that have passed (95780) using the aver- 
age of the best periods from the 2002 and 2003 datasets as 
the mean period. We can then check if we can produce the 
correct phasing using our best periods for the two epochs to- 
gether with the two possible values of /, and their errors. As 
a result we conclude that the data can be phased together 
using either of the two / values. The larger / produces a 
better fit, but we cannot exclude the smaller /. 

The resulting phase folded light curves are plotted in 
Figure 3. The modulation and its asymmetry is roughly 
the same in both datasets. However, it seems like the 2003 
data shows a broader, flatter flux maximum than the earlier 
dataset. The two datasets are identical in terms of instru- 
mentation used. 



5 DISCUSSION 

The current competing models for the system can be classi- 
fied into two types: accreting and non-accreting models (cf 
§1). The accreting and the non-accreting models can be dis- 
tinguished by determining the period change in the system. 

As /// = —P/P, we can calculate / from the orbital 
parameters of the system. For an accreting system, the pe- 
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riod change is given by: 



P 
P 



(Appendix A), where J is the orbital angular momentum, 
M2 is the mass of the secondary star, q is the ratio of the 
mass of the secondary star to the primary star. For systems 
with matter transfer from the secondary to the primary, 
M2 is negative whereas (1 — q) is positive. In general, P/P 
may take a positive or negative value depending on J and 
the driving mechanism of the mass transfer. However, if the 
secondary is a degenerate star, P/P is positive (i.e. negative 
///) when the angular momentum loss from the system is 
due to gravitational radiation (see Ritter 1986). 

For non-accreting system with no coupling between the 
stellar spins and the orbit, 
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P 
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96 / 
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where G is the gravitational constant, c is the speed of light 
M is the total mass of the two stars, and a is the orbital sep- 
aration. If the spin-orbit coupling occurs and the the angu- 
lar momenta of the orbit and stellar spins maintain roughly 
constant ratios, then 



P 



(1 + 9) 2 
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(Appendix A) , where Ri and R2 are the radii of the primary 
and secondary respectively, ot\ and a 2 are the coupling pa- 
rameters of the spins of the two stars with the orbit. For 
parameters typical of double degenerated systems, P/P is 
negative and hence /// is positive (with and without spin- 
orbit coupling) , which is in contrast to the prediction by the 
accreting models. 

If the 321.5 s period is the orbital period of the sys- 
tem, the observed period change is not consistent with the 
prediction of the accreting models but is consistent with 
non-accreting models. If we interpret that the observed pe- 
riod is the spin period of the white dwarf in an intermediate 
polar, both spin up and down is possible. However, we con- 
sider that the spin-interpretation is unlikely (Cropper et al. 
2003). 

In contrast, the detected decrease in the orbital period 
strongly favours the non-accreting models, in particular the 
unipolar- inductor model (Wu et al. 2002), in which the two 
stars in the systems are white dwarfs. In this model, angular 
momentum loss from the system and system evolution is 
due to gravitational radiation and the spin-orbit coupling 
process (if present) via unipolar induction. 

In the abscence of accretion the secondary must reside 
well within its Roche lobe. From the Roche-lobe relation 
of Eggleton (1983) and the Nauenberg (1972) mass-radius 
relation for white dwarfs, we find that for M\ = 1.0 Mq 
and an orbital period of 321.5 s, the white dwarf secondary, 
M2 > 0.14 Mq in order that it stays within its Roche lobe. 

Now, suppose that the observed period is the orbital 
period and that the change is due to orbital evolution. We 
show in Figure 4 the predicted change in frequency for a 
range of component masses, for the case with no spin-orbit 
coupling. Using the longer ROSAT period, we find that for 
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Figure 4. The / as a function of the component masses (assum- 
ing no spin-orbit coupling). The two labelled contours show the 
measured / for the two different ROSAT periods and the contours 
next to them give the 3 sigma limits. 



Mi = 1.0 Mq a secondary mass of 0.3 Mq is required. For 
a less massive primary, a more massive secondary is re- 
quired. Taking the shorter ROSAT period, a Mi=1.0 M 
primary requires a secondary of mass >0.15 Mq. The studies 
of double white dwarf systems (eg Marsh 2000, Napiwotzki 
et al. 2002) have shown that all such systems tend to have a 
mass ratios of the order of 1.0. Furthermore, the individual 
white dwarf masses in these systems are below 0.5 M . Such 
masses are compatible with the measured period change. In 
all of the above mentioned cases, the deduced masses for the 
secondary are consistent with requirement that the star is 
detached from its critical Roche surface. 

If the masses of the system could be determined inde- 
pendently, the degree to which the period decrease differs 
from that expected from gravitational radiation losses alone 
would provide information on the degree of asynchronism 
in the spin of the primary. This asynchronism is the ori- 
gin of the power generation in the unipolar inductor model, 
so continuing timing observations are essential to determine 
whether this is consistent with the observed luminosity and 
reasonable component masses. 

To conclude we note that for perfectly acceptable 
masses of the components, the / is compatible with grav- 
itational radiation. This is consistent with the results for 
V407 Vul (Strohmeyer 2002). The implication is that there 
is no mass transfer in either of these systems, favouring the 
unipolar inductor (electric Star) model. 
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APPENDIX A: ORBITAL EVOLUTION 

Consider a binary with masses Mi and M2 and radii Ri 
and i?2- Let the binary orbital period be P and the orbital 
separation be a. Then u> (= 2tt/P) is orbital angular veloc- 
ity, and M is the total mass. The angular momenta of the 
orbital, star 1 and star 2 are respectively 

, MiM 2 2 

Jo = u, 

2 , 
5 

2 2 

Ji = -M2R2 UJ2 , 

5 

where Wi and L02 are the spin angular velocity of the two 
stars (Assuming a spherical secondary). 

Define two spin parameters ai = uii/u) and Q2 = UJ2/0J. 
The total angular momentum of the binary system can then 
be expressed as 

J = Jo + Ji + J 2 

\ -1/3 1 \ 
= A Ul + A* LO , 

where the two coefficients A and A„ are 
G 2/3 MiM 2 



A 



A. = % (aiMiRl + a 2 M 2 i?I) . 
5 v ' 

If there is no mass exchange between the two stars, dX /dt = 
0. Moreover, if the spin-orbit coupling (via tidal or via 
magnetic field) maintains roughly constant angular momen- 
tum ratios between the orbit and the stars, we may set 
d\*/dt = and obtain a simple relation that relates the 
time derivatives of J and uo: 



J = -- 



,,"/3 



For systems with angular-momentum loss due to grav- 
itational radiation, 



J 



32 G ( Ml Ml 



5 c 5 



M 2 



(see e.g. Landau and Lifshitz 1958), where c is the speed of 
light. It follows that 
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96 / GM\ 3 /c\ 
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The first term above is the ratio of the gravitation radius 
of the binary to the orbital separation; the second term is 
the reciprocal of the light crossing time of the binary; the 
square bracket term can be considered as a structural factor 
of the system; and the last term is the spin-orbit coupling 
correction . For systems with the spins of the stars decoupled 
from the orbit, «i = a 2 = 0, and 
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From above, we can also see that P/P is always negative 
when the {....} term is positive. This can be satisfied easily 
for positive ai and a.2- This implies that the orbit period of a 
binary system decreases because the gravitational radiation 
that it emits carries aways the orbital angular momentum, 
and if in the presence of spin-orbit coupling orbital angular 
momentum is also extracted to synchronise the spins of the 
stars. We note that for systems with zero mass loss but non- 
zero mass exchange between the two stars, 

Jo Mi M 2 la 

X ~ Mi + M2 + 2 a 
M 2 1 w 

= W2 {l - q) -Zu- 



It follows that 

IP 
3P 



Jo M 2 , 
X-M2 (1 - 9) 



(Here, J consists of contributions of orbital angular mo- 
mentum loss by gravitational radiation and by spin-orbital 
coupling.) 



(Mi + M 2 ) 1 /3 ' 



